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Abstract 
(N 

^ \ We present a general analysis that reveals new aspects of the leptogenesis bounds 

Tij- \ on neutrino masses and on the reheat temperature of the Universe. After revisit- 

ing a known effect coming from an unbounded term in the total CP asymmetry, 
we show that an unbounded term in the flavored CP asymmetries has a stronger 



^^ , impact. It relaxes the lower bound on the reheat temperature down to 10® GeV for 

f^ . (M2 — Ml) /Ml = 0(1 — 100) and for a mild tuning of the parameters in the see-saw 

orthogonal matrix. We also consider the effect of the Higgs asymmetry, showing 
\^ . that it lowers the upper bound on the neutrino masses in the so-called fully flavored 

regime where classic Boltzmann equations can be used. Imposing independence of 
the initial conditions contributes to lower the upper bound on neutrino masses as 
well. We study the conditions for the validity of the usual A^i-dominated scenario 
and for the applicability of the lower bound on the lightest right-handed (RH) neu- 
trino mass Ml. We find that except for the two effective RH neutrino scenario, re- 
covered for M3 S> 10^"^ GeV, and for values M2 < 0(10^^ GeV), the final asymmetry 
is more naturally dominated by the contribution from A2-decays. Finally, we con- 
firm in a general way that going beyond the hierarchical limit, the effect of washout 
addition makes the lower bound on Mi more stringent for (M2 — Mi) /Mi = 0(0.1). 



1 Introduction 

Leptogenesis [1] provides an elegant solution to the problem of the non-observation of 
primordial antimatter in the Universe. An appealing feature of this mechanism is that 
it relies on a minimal extension of the Standard Model (SM) where RH neutrinos are 
added to the Lagrangian. In the see-saw limit [2] , small neutrino masses can be naturally 
accommodated in agreement with the data. Considering, for definiteness, the simplest 
and best motivated case of three additional RH neutrinos, only few of the eighteen new 
parameters, three neutrino masses, three mixing angles and three CP violating phases, 
are or can hopefully be probed in neutrino experiments. The remaining 9 'high-energy' 
parameters are out of reach in low energy experiments. In this respect, leptogenesis 
represents also an important cosmological tool to access information on this 'dark side' 
of the see-saw parameter space. 

In an unflavored analysis and in a traditional A'^i-dominated scenario, where it is 
assumed that the final asymmetry is dominantly produced from the lightest RH neutrino 
decays, one finds a lower bound on the lightest RH neutrino mass Mi [3, 4]. This is given 
by Ml > 3 X 10^ GeV [5, 6] at the onset of the strong washout regime, where the final 
asymmetry does not depend on the initial conditions. This lower bound also implies an 
associated lower bound on the reheat temperature, Tj-ch > 1.5 x 10^ GeV [7, 8, 6]. In 
addition an upper bound on the neutrino masses, rrii < 0.1 eV, holds as well [9]. 

A mild hierarchy in the spectrum of the RH neutrino masses, such that M2 > 3 Mi, is 
a necessary condition for the validity of the A'^i-dominated scenario but is not sufficient. 
Indeed, upon closer inspection, the observed asymmetry can still be generated from the 
decays of the next-to-lightest RH neutrinos, realizing a A'^2-dominated scenario [5]. The 
lower bound on Mi is replaced by a lower bound on M2 but this still implies a lower 
bound on T^ch- It is therefore more correct to say that leptogenesis yields a lower bound 
on Trch rather than on Mi. 

Even when flavor effects are considered [10, 11] (see also [12, 13]), the lower bound on 
Treh has been found not to change [6]. On the other hand they induce other important 
modiflcations. First, the flnal asymmetry receives an additional contribution that also 
depends on neutrino mixing parameters [10, 14] and an interesting feature is that such an 
additional contribution can originate solely by low-energy (Dirac or Majorana) phases [10] 
and can even explain the whole observed asymmetry [6, 15, 16, 17, 18, 19]. Second, flavor 
effects can relax the stringent upper bound on the neutrino masses, ttij < 0.1 eV, holding 
in the unflavored regime when a hierarchical heavy neutrino mass spectrum is assumed 
[9]. In [11], it was found that flavor effects completely erase this upper bound. However, 



in [20], it was pointed out that this conclusion rehes on the use of classical Boltzmann 
equations beyond their range of validity and has therefore to be checked within a more 
general description making use of density matrix equations. In [21] it was found that 
using classic Boltzmann equations neutrino masses as large as 2 eV are possible, a value 
larger than the upper bound holding in the unflavored regime and than the current upper 
bound from cosmological observations. A similar conclusion has been recently obtained 
in [22] as well. 

Many of these results have been obtained employing different assumptions or restric- 
tions in the parameter space. In this paper we revisit in detail the leptogenesis bounds still 
assuming classical Boltzmann equations and the A'^i-dominated scenario but without fur- 
ther restrictions on the parameter space and relaxing many assumptions that are usually 
made. For example, we find more general conditions for the validity of the A^i-dominated 
scenario. In general, our analysis reveals various new aspects including effects confirming 
the necessity to go beyond a classical Boltzmann kinetic approach to solve the problem 
of the upper bound on the neutrino masses. 

In Section 2 we set up the general notation. In Section 3 we show that the final 
asymmetry can be written as the sum of different contributions acting independently of 
each other. The first step is to distinguish between a contribution from the lightest RH 
neutrino decays and a contribution from the heavier RH neutrino decays, just two in our 
case. In a traditional A^i-dominated scenario the first one is dominant. This can be still 
conveniently re-cast as the sum of an unfiavored term plus a fiavored term. Finally, the 
unfiavored term can be further decomposed into a piece proportional to a contribution to 
the total CP asymmetry respecting the usual upper bound [23, 3] and into a term that 
is not upper bounded but vanishes when M2 = M3. If M2 7^ M3, this term is typically 
strongly suppressed when a mild hierarchy is assumed and barring a fine-tuned choice of 
the parameters [24, 25, 5]. Analogously, we recast also the fiavored CP asymmetries as the 
sum of an upper-bounded term plus an unbounded extra-term that has been neglected in 
previous works but that proves to be important in interesting cases. 

In the end of Section 3 we review, as a starting point for our analysis, the neutrino 
mass bounds that arise in the minimal scenario, that we call 'vanilla leptogenesis', when 
all simplifications are made at once: the heavier RH neutrino contribution, the fiavored 
terms and the extra-term to the total CP asymmetry are neglected and the hierarchical 
limit is taken. In this case the final asymmetry depends only on 6 parameters [7, 5]. 

In Section 4 we show the role played by the unbounded term in the total CP asymmetry 
of the lightest RH neutrino, Ei, that has to be considered when M2 7^ M3 [24, 25, 5]. 
We study its effect on the bounds in a quite conservative mild hierarchical limit M2 ^ 



3 Ml, showing how a rotation in the complex 2-3 plane is the most effective ingredient 
in enhancing this term and showing that the bounds get modified only for fine tuned 
parameter choices. Our analysis shows moreover that this extra-term acts independently 
of flavor effects. 

In Section 5 we study how the bounds change when flavor effects are taken into account 
and, interestingly, we show that a so far neglected term in the flavored CP asymmetries 
is actually able to relax the lower bound on the reheating temperature in presence of 
wash-out. In the case of the upper bound on the absolute neutrino mass scale, we show 
that this is sensitive to a variation of different assumptions and conditions. We also 
show that, taking into account the Higgs asymmetry then, within the validity of classical 
Boltzmann equations, one cannot say whether the bound mi < 0.1 eV holding in the 
unflavored regime is evaded. This conclusion is strengthened when independence of the 
initial conditions is imposed. 

In Section 6 we consider the contribution from the two heavier RH neutrinos, including 
flavor effects and still assuming the hierarchical limit, showing a (non-trivial) sufficient 
condition for the A^i-dominated scenario to hold. 

Finally, in Section 7, we show how the neutrino mass bounds change going beyond the 
limit of hierarchical RH neutrinos when flavor effects are included. Contrarily to a naive 
expectation, we show that the allowed region, instead of getting enlarged, actually shrinks 
when mild degeneracies for the heavy neutrinos masses are considered. The bounds get 
relaxed only when the heavy neutrino mass degeneracies are much smaller than those of 
the light neutrinos. 

2 General set-up 

Leptogenesis is based on a popular extension of the Standard Model, 
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C = CsM + iNRa^d^Nm - K^^LaNm^ - - MiN^^N^ + h.c. (i = 1, 2, 3, a = e, /x, r), 

(1) 
where three RH neutrinos Njn, with a Majorana mass term M and Yukawa couplings 
h, are added. After spontaneous symmetry breaking, a Dirac mass term m/5 = v h, is 
generated by the vev v = 174 GeV of the Higgs boson. In the see-saw limit, M ^ rriD, the 
spectrum of neutrino mass eigenstates splits in two sets: 3 very heavy neutrinos, Ni,N2 
and iVs respectively with masses Mi < M2 < M3 almost coinciding with the eigenvalues 
of M, and 3 light neutrinos with masses rrii < 1712 < m^, the eigenvalues of the light 



neutrino mass matrix given by the see-saw formula [2], 

m^ = -rriD — m^ . (2) 

Neutrino oscillation experiments measure two neutrino mass-squared differences. For 
normal schemes one has m| — -ml = Aml^^^ and m| — mf = ArrigQi, whereas for in- 
verted schemes one has m| — -ml = Amg^i and m| — mf = Aml^^. For mi ^ ^riatm = 
■\/Am^tm + ^"^soi ~ (0.050 ± 0.001) eV [26] the spectrum is quasi-degenerate, while for 
nil ^ iTT'soi = A/Arrig^j = (0.0088 ± 0.0001) eV [26] it is fully hierarchical (normal or 
inverted). The most stringent upper bound on the absolute neutrino mass scale comes 
from cosmological observations. Recently, quite a conservative upper bound, 

mi<0.2eV (95% CL) , (3) 

has been obtained by the WMAP collaboration combining CMB, baryon acoustic oscilla- 
tions and supernovae type la observations [27]. 

With leptogenesis, this simple extension of the Standard Model is also able to explain 
the observed baryon asymmetry of the Universe [27] 

r^^^^ = (6.2 ± 0.15) X 10-1° . (4) 

It is widely known that, in order to generate a baryon asymmetry in the early Universe, 
one needs to satisfy the three Sakharov conditions [28]. At temperatures T > 100 GeV, 
baryon number is violated by the non-perturbative sphaleron processes [29]. Moreover 
CP is violated in the decays of the heavy RH neutrinos. Indeed the Dirac mass matrix 
is in general complex and this provides a natural source of CP violation. At the same 
time departure from thermal equilibrium occurs in the decays of the heavy RH neutrinos 
as well. This can be conveniently quantified in terms of the decay parameters, defined as 
Ki = Ti/HT=Mi) given by the ratio of the decay widths to the expansion rate when the 
RH neutrinos become non-relativistic. The decay parameters can be expressed in terms 
of the Yukawa couplings by 

Ki = — , where m^ = — (5) 

m^ Mi 

are the effective neutrino masses and m^ is the equilibrium neutrino mass [7] given by 
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m^ = ^^-^ ~ 1.08 X 10-^eV. (6) 

3^5 Mpi ^ ' 



There are two ways how CP violation can manifest itself. A first one is given by having a 
decay rate of Ni into leptons, Fj, different from the decay rate into anti-leptons, Fi. This 
is parameterized by the total CP asymmetries 

F - F 

A perturbative calculation from the interference of tree level with one loop self-energy 
and vertex diagrams gives [30] 



167r ^ {Wh),, v^^^T^ 

having introduced [9] 

■l + a;\ 2-x 



i{x) = -x 



x) In 



X I I — X 



(9) 



A second way [10] is given by the possibility that in general, indicating with |£j) the 
final lepton quantum state and with |£^) the final anti-lepton quantum state, one has 
l^j) 7^ CP\C'^. This can be easily understood when the flavor composition of the final 
lepton state is considered. Indeed, introducing the projectors on the flavor eigenstates, 
they can be written like the sum of two terms, 

p.« = m\Q? = pl + ^ (10) 

P^. - m\Q? = PL-^- (11) 

The first term is the tree level contribution and is common to both projectors but the 
second term, from loop corrections, changes sign and gives rise to a different flavor com- 
position when APja 7^ 0. These two CP violating effects translate respectively in two 
separate terms in the flavored CP asymmetries, 

F — F A P 

_ ^ la i- ta pO I ^-^ ta /-|r,\ 



Ti + Ti 

where Fjq, = Pjq,Fj and Fj,^ = Pia^i. The flavored CP asymmetries can be calculated 
using [30] 

(13) 



The second CP violating contribution yields an additional contribution to the final asym- 
metry only if the fiavor composition plays a role in the determination of the final asym- 
metry. This depends on the effectiveness of the charged lepton interactions implied by the 
term fa^La^Ra^ in the Lagrangian, which is diagonal in flavor space. The latter implies 
that the processes ia(^a ^^ ^ and ^qCq, ^^ $ A^ and the CP conjugated, (where A^ (a = 
1, 2, 3) are the SU{2)l gauge bosons) occur at a rate Fq ~ 5 x 10"^ T f^{a = e, fj,, r) [31]. 
If these processes are effective, then they measure the flavor composition of the flnal lep- 
tons and this becomes a relevant ingredient in the determination of the final asymmetry 
if a second condition is fulfilled as well, as we will comment. 

In the hierarchical limit the decays of just one species of RH neutrino, the lightest 
[1] or the next-to-lightest [5], dominantly contribute to the final asymmetry. If Fq, <^ 
FJj-, {i = 1,2) during the relevant period of the asymmetry generation, where F^p denotes 
the inverse-decay rate, then the coherence of the lepton states is preserved on average 
between a decay and a subsequent inverse decays and the unfiavored regime, where fiavor 
effects are negligible, holds. This requirement implies [20] 

Mi>5xlO"GeV. (14) 

In this case an approximate set of Boltzmann equations is given by 



dz 

dNB-L 

dz 



-A (iViv, - iV^^) , 2 = 1,2,3 (15) 

3 

Y^ e. A {Nm. - N2) - Nb-l [^W{z) + J] W^z)] , (16) 



=1 



where z = Mi/T and where we indicated with Nx any particle number or asymmetry 
X calculated in a portion of co-moving volume containing one heavy neutrino in ultra- 
relativistic thermal equilibrium, so that N'^,{T ^ Mj) = 1. With this convention, the 
predicted baryon-to-photon ratio rjB is related to the final value of the final B — L asym- 
metry by the relation 

VB = «sph-^ ^ 0.96 X 10-2iV|,_^ , (17) 
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where A^^*"^ ~ 37, and Cgph = 28/79. Defining Xi = Mf/Mf and Zi = z i/a^, the decay 
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factors are given by 

n — -'^D.i r^ /I 

Hz \7, 



where H is the expansion rate. The total decay rates, F^ j = Fj -|- Fj, are the product of 
the decay widths times the thermally averaged dilation factors (I/7), given by the ratio 



JCi{z) / IC2{z) of the modified Bessel functions. Tlie equilibrium abundance and its rate 
are also expressed through the modified Bessel functions, 

After proper subtraction of the resonant contribution from AL = 2 processes [32], the 
inverse decay washout terms are simply given by 

Wl''{z) = ^K,y^,ICi{z,)zl (20) 

The washout term AW{z) is the non-resonant AL = 2 processes contribution. It gives a 
non-negligible effect only at z ^ 1 and in this case it can be approximated as [7] 

where u ^ 0.186 and m^ = ml + m^ + m\. Notice that we are neglecting AL = 1 
scatterings [33, 14], giving a correction to a level less than ~ 10% [6], thermal corrections 
[8], giving relevant (though with big theoretical uncertainties) corrections only in the weak 
washout, and spectator processes [34, 35], that produce corrections to a level less than 
~ 20% [35]. 

On the other hand, if the charged lepton Yukawa interactions are in equilibrium (Fq, > 
H) and faster than inverse decays, i.e. 

r« > rJD , (22) 

during the relevant period of the asymmetry generation, then lepton quantum states lose 
coherence between the production at decay and the subsequent absorption in inverse 
processes. If the quantum state becomes completely incoherent and is fully projected 
on one of the flavor eigenstates, each lepton flavor eigenstate £q, can be treated as a 
statistically independent particle species and a 'fully flavored regime' is obtained. Note 
that one has to distinguish a two-flavor regime, for Mj > 10^ GeV, such that the condition 
Eq. (22) is satisfied only for a = r, and a three-flavor regime, where it applies also to 
q; = /i. 

In the fully flavored regime (two or three flavors), classical Boltzmann equations can be 
still used like in the unflavored regime, with the difference, in general, that now each single 
flavor asymmetry has to be tracked independently. Since sphaleron processes conserve the 
quantities Aq, = 5/3 — Lq, (a = e,/i, r), these are the convenient independent variables 



to be used in the set of Boltzinann equations that can be written as 

dNN^ 



dz 
dz 



-A (iV^. - iV^^) (^= 1,2,3) (23) 

Y. e.^ A {N^. - KV - E ^^ (^^/3 + C^) W^/° N^, , (24) 



where we are using the same approximations as in the unflavored case but neglecting the 
non-resonant AL = 2 term, since this counts only for Mi > 10^^ GeV (^^atm/ Yli ^^f ) ^^^^ 
also a contribution from AL = processes that one has to consider in the flavored case. 
Notice that the final B — L asymmetry is now calculated as N^b-l — Yla ^a^ ■ 

The C^ matrix [12] relates the asymmetries stored in lepton doublets, i^, to the asym- 
metries Afj = B/3 — La asymmetries and is given, in a two-fiavor regime, by [14] 

C'=^(^^^ -^^] (25) 

316 \^ -17 208 J' ^ ^ 

The C^ matrix takes into account the washout due to the asymmetry stored in the Higgs 
field [10] and is given by 

C^ = ^(41,56). (26) 

The Higgs asymmetry has been neglected so far but, as we will see, it has a relevant effect 
on the upper bound on the neutrino masses. Indeed the sum of the two matrices gives 

C.C' + C.f"! °^2^V (27) 

y 0.21 1.01 J 

The off-diagonal terms give a small effect in the calculation of the final asymmetry [14] 
and therefore, in the end, one can safely use the approximation C ~ / in the derivation 
of the bounds. In Section 5 we will compare the results when the Higgs asymmetry is 
neglected, when it is taken into account and when one uses the approximation C ~ /, 
showing that the latter works very well justifying its use for the remainder of the paper. 
Taking for simplicity the two fiavor case, it is instructive to sum over the fiavor Eq. (24), 
obtaining 



dN 



r-^ 



^ ., Z^ (iV^, - iV^^) - i iV^_, X^ W^i° + i [iV^„ - AT^ J ^ (i^^ - i^y W^^ 



dz ^ ' ' ''' ''^' 2 ^"" ^ ' 2 

1=1 i=l 4=1 

(28) 
This equation clearly shows that when the washout vanishes there must be no difference 
between the unflavored and the fully flavored regime. Therefore, the lower bounds on 
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Ml and on T^^h obtained in the limit of no-washout, rigorously for an initial thermal 
abundance and approximately for an initial vanishing abundance, do not change when 
flavor effects are taken into account [6]. 

A convenient parametrization of the Dirac mass matrix is obtained in terms of the 
orthogonal matrix [36] 

mD = UDy,^nD'^\ (29) 

where we defined Dm = dia,g{mi,m2.,m^) and Dm = diag(Mi, M2, M3). The matrix U 
diagonalizes the light neutrino mass matrix rriy, such that U"^ iriyU*^ = —Dm, and it can 
be identified with the lepton mixing matrix in a basis where the charged lepton mass 
matrix is diagonal. We will adopt the parametrization [37] 



U 



( C12C13 S12C13 sise '^ \ 

-S12C23 -Ci2S23Sl3e'^ C12 C23 - S12 S23 Sl3 e"^ S23C13 



. sf]^ . '±'2 

X diag(e'~,e'~, 1) , 

\ S12 S23 ~ C12 C23 Si3 e' — C12 S23 "~ -512 C23 Si3 e' C23 C13 j 

(30) 
where Sij = sin^jj, Cjj = cos 6ij and, neglecting the statistical errors, we will use 612 = 
7r/5 and ^23 = 7r/4, compatible with the results from neutrino oscillation experiments. 
Moreover, we will adopt the 3o" range S13 = — 0.20, allowed from a global Su analysis for 
unitary U [26] , an approximation that holds with great precision in the see-saw limit with 
Mi ^ 100 GeV. With the adopted convention for the light neutrino masses, mi < 1712 < 
m^, this parametrization is valid only for normal hierarchy, while for inverted hierarchy 
one has to perform a column cyclic permutation. In a general analysis, leptogenesis 
bounds are not depending on the scheme, normal or inverted, but in restricted scenarios, 
like in the effective two RH neutrino scenario where the third is very heavy and decouples 
or in 'Dirac phase leptogenesis' [19], differences can arise and depend on flavor effects. 
We will signal these differences in our analysis. 

It will also prove useful to introduce the following parametrization for the see-saw 
orthogonal matrix in terms of complex rotations 



^(c^2i, c^3i, ^^32) = diag(±, ±, ±) i?i2(c^2i) ^13(^31) ^23(^32) , (31) 



where 



R12 ■ 



v/T^^ -C..1 ^ / iyT^T^ -^31 ^ / 1 « ^ 



V 01 



1 \ , R 



±\/l - t^i2 -"32 



V "31 ±\/l-"|i / \ CJ32 ±^/T^i^'i 

(32) 
and where the overall sign takes into account the possibility of a parity transformation as 
well. Notice that, using the orthogonal parametrization, Eq. (29), the effective neutrino 
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masses, and consequently the decay parameters, can be expressed as linear combinations 
of the neutrino masses [38, 7], such that rhi = ^ . rrij |i^^j|. Notice that the orthogonality 
of Q is equivalent to the see-saw relation for the light neutrino masses and in particular 
one has that Re[r2fj] is the contribution to rrii from the term oc 1/Mj. Therefore, large 
absolute values of the Q entries imply a strong fine tuning not only because they require 
phase cancelations but also because they imply that neutrino masses are much lighter than 
terms oc m|,/M because of sign cancelations. Therefore, such choices tend to transfer the 
explanation of neutrino lightness from the see-saw mechanism to some other mechanism 
that has to explain the fine-tuned cancelations. The interest for considering models with 
very large \^ij\ is merely phenomenological since they make possible to satisfy neutrino 
masses with the see-saw mechanism and at the same time to have TeV RH neutrinos 
with large Yukawa's, making possible to detect them in colliders [39]. This will not be 
our point of view in this paper and we will conventionally consider orthogonal matrices 
to be 'reasonable' if \ujij\ < 1, implying \^ij\ < 1, and 'acceptable' if \ujij\ < 10, implying 
\flij\ < 10. These will be the two benchmark cases that we will adopt in the plots. 

To conclude this Section we just want to recall the particularly relevant case of two 
effective RH neutrinos, obtained in the limit where M3 ^ 10^^ GeV [40, 41, 42]. In this 
limit the orthogonal matrix necessarily collapses into 

/ 1 \ 

n= ±v/l-fiii -^31 , (33) 

corresponding to have UJ32 = 1 and UJ21 = 1 in the Eq. (32). 

3 Vanilla leptogenesis 

Using the approximation C c::^ I, & general solution for the final asymmetry can be written 

as [6] 



I, a 



with the final values of the 9 efficiency factors given by 



^ (K- P° ) = - r dz' ^^^' e-i^w{z)+j:^ Pl /J dz"wl^(z",m ^ 



kUK,, Pi) = - dz' ^ e-l^^^^^+^» ^- J^' '^^ ^r(-";^.)J . (35) 

This solution holds both in the fully flavored regime and in the unfiavored regime, adopting 
the convention that when the condition (14) applies, all projectors Pj^^ = 1. It is indeed 
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easy to verify that in this case, summing over the flavor, the set of equations (24) reduces 
to Eq. (16). On the other hand it should be also noticed that Eq. (24) holds only if the 
condition (22) is respected and this, when applied to the A^^i decays and inverse decays, 
translates into 



Ml < ttT77!d7TT7^^- (36) 



10^2 GeV 

For WI^{zb) > 1 there is an intermediate regime where the unflavored regime does not 
hold but at the same time the validity of the Eq. (24), and therefore of the expression (34), 
is not guaranteed. We will signal in the plots the results obtained in the fully flavored 
regime but for which the condition (36) is not satisfled. These results should be therefore 
checked within a more general kinetic description employing density matrix equations 
able to describe the regime where coherence (or decoherence) of the final lepton quantum 
state is only partial. 

Notice that the N^^ 's are the values of possible pre-existing flavored asymmetries. 
In the unflavored regime Ki is the only parameter that determines whether the flnal 
asymmetry depends or not on a possible pre-existing asymmetry [9] . Taking into account 
flavor effects the problem is more involved and there are different issues to be considered. 
Here we will not face this problem and we will simply assume that the flrst term is 
negligible. 

It is instructive to re-cast the Eq. (34) in an approximate way that enlightens different 
effects and contributions. For deflniteness, we consider the two-flavor regime holding for 
Ml > 10^ GeV. The total CP asymmetry ei can be recast as [5] 



^1 = 


= ^{x2) ei{Mi, mi,uj2i, ^31) + [^(a^a) - ^(2:2)] A£i(Mi, mi, Q) , 


(37) 


where, deflning 






e(Mi) E 


- T^ 2 ^^^ (3{mi,U2i,iU3i) = ^^F^ p^^, 

167r v^ matm l^j rrij iVj^\ 


(38) 


one has 








£i(Mi, mi, uj2i,uj3i) = e{Mi) /5(mi, ^21, ^31) 


(39) 


and 


A.i(Mi,mi,l])-e(Mi)'"^[^'^"^^?^^^^]\ 


(40) 



matm mi 
With these deflnitions, the flnal asymmetry can be written approximately as 

Nl,_L ~ N{i{^{x2)eiiMi,rm,U2i,iU3i) + Ux^)-^ix2)]AeiiMi,rm,n)} k{ (41) 
, APi^(Mi,mi,(],^) ^ 1 fi^Y-r f^ fi 

a 
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where Nf[ is an effective number of ffavors. In the unffavored regime one has Nf[ = 1, and 
in the fully flavored regime, if both flavors experience a strong washout (P^^^-ft'i ^ 1), 
one has approximately Nq ~ 2, while in general 1 < A^h ^ 2. 

The simplest scenario, that we call vanilla leptogenesis borrowing the name from ob- 
servational cosmology, corresponds to taking all possible simplifying assumptions: 

1. hierarchical limit, M2 > 3 Mi, so that ^{x2) — 1; 

2. negligible contribution from the heavier RH neutrinos; 

3. negligible flavor effects (A^a = 1 and APia = 0); 

4. M2 = M3, so that ^{x2) - ^(xg) = 0. 

The calculation of the final asymmetry then simply reduces to 

Nb-l ^ ei{mi, Mi,uj2i,uj3i) K{{mi, Mi,uj2i,uj3i) , (42) 

depending only on 6 unknown parameters. The efficiency factor is well approximated 

by [7] 



K{{M,,mi,K,) = k{{K,) exp I -^ (- 



Ml \ / m ^ ^ 



IQio GeVy VeV 
For the case of a thermal initial A^i-abundance (A^j^^ = 1), one has 



^{{Ki) ~ k{Ki) = ^ 



1-exp ( --KiZb{Ki] 



(43) 



(44) 



Kiz^iKi 
where zb is approximately given by the expression [43] 

ZBiKi) ^ 2 + 4iri°-i=' exp (-|^) • (45) 

In the relevant range 5 < -ft'i < 100 this expression is further well approximated by the 
simple power-law ^^(i^i) ~ 0.5/ K}'^ [44]. 

In the case of a vanishing initial A^^i -abundance {N^^ = 0), one has to take into account 
both a negative and a positive contribution, such that 

k{{Ki) = Ki{Ki) + k'^Ki) . (46) 

The analytic expressions for k,^_{Ki) and k,^_^{Ki) can be found in [7]. Imposing that 

the predicted final asymmetry, Eq. (17), explains the observed one, Eq. (4), yields the 

condition 

M 

K\{Mi, mi, Ki) fJ{mi, u;2i, u;3i) 
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Figure 1: Vanilla leptogenesis and leptogenesis conspiracy. The allowed region in the 
plane (mi, Mi) is shown for different values of rriatm setting the value of M (cf. Eq. (48)). 
In the central panel the true measured value is considered, wiatm — 0.050 eV, while in 
the left and right panels two 'wrong' values are considered, matm = 10~^eV (left) and 
^^atm = 10 eV (right). In all three panels, the red crosses correspond to a thermal initial 
iVi-abundance while the green ones to a vanishing initial A'^i-abundance. The hatched 
area indicates values of mi excluded by the upper bound Eq. (3). 



where we defined 



— 167riV^ 



rec 2 ^CMB 
^ 'IB 



'^sph '''■atm 



(6.6 ± 0.3) X 10*^ GeV > 5.7 x 10*^GeV 



(4J 



The last inequality gives the 3a value of M that we used in the plots to obtain the allowed 
region in the (mi. Mi) plane scanning over all values of ti;2i and u^i. The result is shown 
in the central panel of Fig. 1. One can notice the presence of the usual lower bound on 
the lightest RH neutrino mass. Mi > 2.3 x 10^ GeV for initial vanishing abundance and 
Ml ^ 5.7 X 10^ GeV for initial thermal abundance [4]. This can be easily inferred from 
the Eq. (47) neglecting the exponential factor in Eq. (43) and using the well-known upper 
bound 



P{mi,U2l,UJ3i) < 



'^atm 



f {nil, mi) 



(49) 



mi + 1713 

where the function < f {1711,7711) < 1 is unity in the limit fhi/mi -^ 00, and vanishes 
for mi = rrii. This function can be derived analytically together with simple analytic 
expressions valid in particular regimes [5]. For mi <^ mi <^ m^ and for initial thermal 
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Figure 2: Allowed region in the plane {mi, Mi) for fi = i?i3 (left panel), and Mi lower 
bound versus Ki for the effective 2 RH neutrino case obtained for M3 — i> 00 and implying 
mi = (right panel). In the left panel, the color coding is the same as in the previous 
figure. In the right panel, the red crosses correspond to normal hierarchy whereas the 
green ones to inverted scheme. 

abundance one simply finds Mi > M > 5.7 x 10^ GeV, in agreement with the numerical 
result. From the central panel of Fig. 1 one can also notice that vanilla letogenesis predicts 
mi < 0.12 eV [9], in agreement with the current observational upper bound (cf. Eq. (3)). 
In the plots the hatched region indicates the excluded values. This upper bound can be 
derived analytically as well [7]. 

We performed a simple exercise showing how the allowed region in the {mi , Mi ) plane 
would have been for values of rnatm and m,soi different from the true measured ones. We 
kept the ratio matm/'^soi constant. In the left panel matm = 10~^ eV, while in the right 
panel rriatm = 10 eV. One can see how the allowed region shrinks considerably, and 
almost disappears for these extreme values. This is one way to show the 'leptogenesis 
conspiracy' [45], that means how, order-of-magnitude-wise, the measured atmospheric 
and solar neutrino mass scales are optimal for leptogenesis to be successful. It should 
be noticed that for ?7iatm = 10 eV, even though the lower bound on Mi is much lower, 
the density of points in the allowed region is very low since they correspond to a very 
fine-tuned situation where fhi ^ matm- 

In Fig. 2 we show the bounds for particular choices of the orthogonal matrix: in the 
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left panel for Q = R13, corresponding to CJ21 = ti;32 = 0, and in the right panel for the 
effective two RH neutrino case obtained in the limit M^/IO^'^ GeV ^ 1 and corresponding 
to an orthogonal matrix with a;2i = ^32 = 1. In the first case, one can see how the bounds 
do not change compared to the general case, showing that this is the choice saturating the 
bounds, a well-known result [4] . Notice that for Q = R13 there is no dependence of the 
bounds on the light neutrino mass scheme, normal or inverted, and therefore, in vanilla 
leptogenesis, normal or inverted schemes produce the same bounds [4]. On the other 
hand, in the second case the lower bound on Mi becomes more stringent [41], especially 
in the case of inverted hierarchy. Indeed, in the vanilla case, for the same choice of the 
orthogonal matrix, the final asymmetry in an inverted scheme can be only less than in a 
normal scheme, or at most equal in the special case Q = -R13, as discussed analytically 
in [5]. This result is easy to understand qualitatively: the dominant term in the CP 
asymmetry is suppressed when the neutrino masses increase, either when mi increases, 
or switching from normal to inverted hierarchy since in this case m2 gets higher. 

In the next sections we will relax the assumptions of vanilla leptogenesis, studying how 
the leptogenesis bounds change accordingly. In most cases the effect of the assumptions 
on the bounds is independent of each other and therefore they can be studied individually. 
However, in a few cases the interplay of different effects can yield interesting interferences. 
For example, going beyond the hierarchical limit one has also to consider the contribution 
of the heavier RH neutrinos to the final asymmetry and to the washout. 

4 Extra-term in the total CP asymmetry 

In this Section we relax the assumption M2 = M3 defining vanilla leptogenesis, studying 
the effect on the bounds of the term proportional to Aei{Mi,mi,Q) in the Eq. (41), 
and comparing our results with those obtained in [24]. This effect clearly saturates to a 
maximum when M3/M2 ^ 1 and in the plots we fixed M3/M2 = 100. 

From the Eq. (41) it can be noticed that this effect acts independently of flavor effects 
and can be even dominant when the parameters are properly tuned. In particular this 
extra-term is able, although with quite a strong fine-tuning, to relax the lower bound on 
Ml. A separate analysis is fully justified since there is no interference between the two 
effects. In the 6 panels of Fig. 3, we show the crucial role played by the parameter 0032 
[5] in enhancing the extra-term in the total CP asymmetry (cf. Eq. (40)). This relaxes 
the mass bounds in a remarkable way if |a;32| ^ 1 since the CP asymmetry enhancement 
is not counterbalanced by an increase of the washout that is driven by Ki and that is 
independent of CJ32. It can be seen in the top right panel that when 0032 = the bounds 
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Figure 3: Effect of the extra term [cf. Eqs. (40) and (41)] in ei on the neutrino mass 
bounds. In all panels M3 = 100 M2 and M2 = 3 Mi, except in the top center panel where 
M2 = 10 Ml. Top left and top center panel: all three \ujij\ < 1; top right panel: |co'32| = 
while |ti;2i|, \uJ3i\ < 10; bottom left panel: |ti;2i|, 1^7321 < 1, {uj^il < 10; bottom center panel: 
If^sil) 1^32! < I5 |^2i| < 10; bottom right panel: |u;2i|, \uJ3i\ < 1, 1^732! < 10. 
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are almost unchanged, even though we allowed |a;2i|, |u;3i| < 10. 

Notice that Aei increases with the neutrino masses, contrarily to £i. This is the 
reason why it tends to relax the upper bound on mi. It also tends to be higher for 
inverted schemes compared to normal schemes, even though the bounds are saturated 
for a choice of the parameters where there is no dependence on 7712 so that inverted and 
normal schemes give the same results. 

The results of the panels can be easily understood analytically using Eqs. (40) and (41). 
For example, since $,{x) ~ 1 + 5/{9x) when a; ^ 1, the extra-term is suppressed like 
{M1/M2Y. It should be also noticed that the extra-term vanishes exactly in the limit of 
two effective RH neutrinos, obtained for Ma/lO^^ GeV — > 00. 

In conclusion, the possibility to exploit the extra-term oc A£i(Mi,'mi, ^2) to relax the 
lower bound on Mi relies on models where \uj^2\ ^ 0.2 (M2/Mi)^. Therefore, already for 
M2 > 3 Ml, quite a high level of fine tuning is required. Concerning the upper bound 
on iTii the conditions are less stringent: \oJ^2\ ^ 0.2 {M2/M1Y and/or |u;2i| ^ {M2/M1Y, 
such that one has to impose, more conservatively, M2 > 10 Mi. We will see in the next 
Section that flavor effects have a bigger impact in relaxing the bounds compared to the 
vanilla scenario. 

5 Adding flavor to vaniUa 

Relaxing only the second assumption defining vanilla leptogenesis, one has 



iVfl£i(Mi,mi,cu2i,cc^3i) 4(Mi,mi,7^i) + - APi^{Mi,mi,Q,U) [k{^ - ^y . 



1 
2 

From Eq. (13), in the HL, one finds [6] Eia = ^la + ^-^^la, where 

3 



^la 



i6.(ftt„,);Tg ;!=""[''"'"-<"'">"] <''> 



and 

Taking advantage of the orthogonal parametrization (cf. Eq. (29)) and defining ri^, = 
eia/e{Mi) = fia + Aria, one has 



ria = - V '^i^'"''"" lm[Uah Ki ^hi ^11] (53) 

^ mi matm 
h,l 



and 






' hlk " -^ ^^""^ ^^atm 



The second term has been neglected in previous analyses but, as we will see, it can 
dominate under some conditions relaxing the leptogenesis bounds holding in the vanilla 
scenario. The most important difference between the two terms is that the first is upper 
bounded [11], 

ric< yPio^ms/m^i^, (55) 

while the second is not. From this point of view the A^iq, term is analogous to the extra 
term in the total CP asymmetry but, as we will see, it affects the bounds in a more 
relevant way. 

The analytical expressions for the k\^, generalizing those for k{, can be found in [6]. 
The generalization of the expression (47) becomes 

M 

In the fully flavored regime this gives rise to a lower bound on Mi that always falls in the 
two-flavor regime with negligible washout from AL = 2 processes. Therefore, this can be 
expressed like 

Ml = ; ; ~ J -, J J . (57) 

ri.«:L + ri,e+M4e+M ^fl 4(^i) + | [APi^MMi)] [^^ - 4, J 

From the approximate expression, one can see once more that the lower bound, can be 
relaxed compared to the unfiavored case only if there is some washout, otherwise A^a = 1 
and K,{^ — i^i^e+n = 0- This implies Ki > 1. In the limit of no washout, for Ki <^ 1, one 
recovers the usual lower bound Mi > M in the case of initial thermal abundance. In the 
strong washout a big relaxation is possible only in the one-flavor dominance case, where 
one of the projectors P°^ ^ P^^, otherwise close to the democratic case, where P^^ ~ P^j^, 



the difference [k,\^ — k,\ g , ] tends to suppress the final asymmetry. 



5.1 Lower bound on Mi 

In Fig. 4 we show the allowed region in the (-ft'i,Mi) plane for rrii = 0. The plots are 
obtained scanning the seven free parameters. Mi and the 6 parameters in Q, showing 
only the points where rjB > 77^^^ at 3 a. This is equivalent to search for the points where 
Ml is larger than the right-hand side of the Eq. (56) with M ^ 5.7 x 10® GeV. In this 
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Figure 4: Lower bound on Mi versus Ki for rrii = and imposing \u}ij\ < 1. We show 
the results obtained in the fully flavored regime (red points) comparing them with those 
obtained when they are neglected (green points). Top- left panel: thermal initial A^^i- 
abundance (iV]^^ = 1). Top-center panel: vanishing initial A^i-abundance (iV]^^ = 0). 
Top-right panel: points falling in the strong washout regime where the final asymmetry 
depends on the initial A'^i-abundance to a level less than 10%. All bottom panels assume 
a thermal initial A'^i-abundance. Bottom-left panel: like top-left but removing all points 
that violate the condition (36). Bottom-center and right panel: inverted hierarchy keeping 
or removing points violating the condition (36) respectively. 
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subsection we are imposing \uJij\ < 1, implying an upper bound on Ki. We will study in 
Section 5.3 the effects of turning on large values of \oOij\. 

In the plots the red region is the additional part of the allowed region due to flavor 
effects within the fully flavored regime while the green region is what one obtains within 
vanilla leptogenesis neglecting flavor effects . In the top-left panel the flnal asymmetry 
has been calculated for an initial thermal A'^i-abundance, while in the top-middle panel 
for an initial vanishing A'^i-abundance. One can see again that flavor effects can relax the 
lower bound only in the presence of washout, that means when fi'i > 1 and the amount of 
the relaxation increases with Ki. Essentially the lower bound we flnd coincides with the 
lower bound found analytically in [6] that corresponds to neglect Ari^ and maximizing 
fia and k\^ in the Eq. (56) in the case of one flavor dominance. In [6] the lower bound was 
numerically calculated only for a particular case, Vt = R13, and just a small relaxation 
was found for vanishing rrii. Here, allowing for a;2i 7^ 0, we flnd a large relaxation 
also for vanishing mi. In the top-right panel, we selected only the points for which 
there is independence of the initial conditions, more exactly those for which there is a 
difference to a level less than 10% between thermal and vanishing initial A^i-abundance. 
This generalizes the deflnition of strong washout regime when flavor effects are taken into 
account. The critical value of Ki increases from ~ 3 in the unflavored case to ~ 7 in the 
flavored case. 

In the bottom left panel we flnally show only the points that further satisfy the condi- 
tion of validity of the classic kinetic equations in the fully flavored regime (cf. Eq. (36)), 
and one can see how there are no differences between the unflavored and the fully flavored 
regime since these arise for larger Ki values. 

In the bottom center and right panels we show the results for inverted hierarchy 
without imposing the condition of validity and imposing it, respectively. One can see 
that the situation is not very different from the normal hierarchy case. Slightly larger 
values of Ki are allowed, up to about 2 x i^atm — 93. This implies that the relaxation of 
the lower bound on Mi can be as large as one order of magnitude for the maximal allowed 
value of Ki. Notice that these points still satisfy the condition of validity (36). 

In Fig. 5 we consider two special cases for the orthogonal matrix when mi = 0. In 
the left panel Q = R13 (i.e. a;2i = ^32 = 0). In this case it is easy to show [6] that 
APia = and therefore in the Eq. (50) only the flrst term survives and all flavor effects 
reduce to an enhancement of the flnal asymmetry compared to vanilla leptogenesis given 
just by N{ ~ 2. In the right panel we consider the case of very large M3 ^ 10^^ GeV to 
be compared with the right panel of Fig. 2 in the unflavored case. One can see that this 
time inverted hierarchy (green and purple areas) is not so suppressed compared to normal 
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Figure 5: Lower bound on Mi versus Ki for Vt = R13 (left) and for M3 ^ 10^^ GeV 
(right). In the right panel the red and blue regions correspond to normal hierarchy and 
the purple and the green to inverted hierarchy. The blue and purple regions are obtained 
switching off the PMNS phases. 

hierarchy (red and blue areas) as in the unflavored case. This is due to the presence of 
the PMNS phases that give a further contribution to the asymmetry. Indeed when phases 
are switched off (purple and blue areas) again the allowed region for inverted hierarchy is 
strongly reduced compared to normal hierarchy, similarly, though to a minor extent, to 
what happened in the unflavored case. In this specific case M3 ^ 10^^ GeV, the effect of 
phases has been recently studied in [46]. The fact that PMNS phases can give a dominant 
contribution to the final asymmetry was first noticed in the case Q = R13 and mi 7^ in 
[6]. 

5.2 Upper bound on mi 

In this subsection we allow rrii 7^ 0, investigating how the upper bound on mi and its 
dependence on Mi changes when flavor effects are included. For this purpose, we show 
plots in the {Mi, mi) plane. As in the last subsection, we impose the condition \u}ij\ < 1. 
In all figures we distinguish three different kinds of points: those characterized by 
a strong one- flavor dominance for which P^^ < 0.1 or P°g + P? < 0.1 (red crosses), 
those characterized by a mild one-flavor dominance for which 0.1 < P°^ < 0.45 or 0.1 < 
Pi^ + Pi < 0.45 (green x) and finally those for which a democratic scenario is realized 
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Figure 6: Mi vs. mi for normal (left panel) and inverted (right panel) hierarchy. The 
red straight crosses denote a projector P^^ < 0.1 or P°g + P^^ < 0.1 , the green x's 
0.1 < P^^ < 0.45 or 0.1 < Pi°, + P^^ < 0.45 and the blue stars 0.45 < P^^ < 0.5 or 
0.45 < P°g + P°, < 0.5. 



such that 0.45 < P°^, Pi'e + Pi)j < 0.55 (blue stars). This will make possible to understand 
under which circumstances the upper bound on mi can be evaded when flavor effects are 
taken into account. 

In Fig. 6 the results are shown both for normal (left panel) and inverted (right panel) 
hierarchy and have been obtained using the approximation C = L One can see that there 
is no upper bound on mi, as first pointed out in [11]. The evasion of the bounds occurs 
in a one-flavor dominance, as expected. Note that the results have been obtained without 
imposing the condition of validity of the fully flavored regime Eq. (36). In the case of 
inverted hierarchy (right panel) the bounds do not change significantly. 

In Fig. 7 we consider the special case Q = P13, corresponding to UJ21 = UJ32 = in the 
Eq. (32). The number of free parameters gets therefore reduced to 6 (one of the PMNS 
phases, $2? is irrelevant in this model). One can see that the bounds are very similar 
to the general case, just slightly more stringent at large mi. Therefore, the special case 
Q = Pi3 gives an approximate condition for the saturation of the bounds in the {m,i, Mi) 
plane. The phases in the PMNS matrix can play a crucial role in the fully flavored regime. 
This was first emphasized in [10] and then analyzed in more detail in [14, 6]. Here we 
want to show the effects of the PMNS phases on the lower bound found in the previous 
figures comparing the previous results with those obtained when the PMNS phases are 
turned off. The result for a general Q matrix is not shown because it is given precisely 
by Fig. 6. There are indeed enough phases present in the fl matrix in this case to realize 
a strong one-flavor dominance and to saturate the general lower bound even though the 
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Figure 7: Same as in Fig. 6 but for the special case Q = Ri^. 





Figure 8: Same as Fig. 6, for Q = Ri^ with all PMNS phases turned off. 



PMNS phases are set to zero. The situation is different when one considers special cases 
like Q = Ri3 as first found in [6]. The result is shown in Fig. 8 and comparing with 
Fig. 7 one can see clearly that the PMNS phases are responsible for the saturation of 
the bounds at large mi. The reason is that there is in this case only one phase in the 
Q matrix, and this is not enough to fulfill the conditions for the saturation of the lower 
bound. Therefore, the Majorana phase $i and the Dirac phase 6 play a crucial role in 
this case analogously to the case M3 ^ 10^^ GeV for inverted hierarchy, as seen in the 
previous subsection and recently pointed out in [46]. 

Let us now discuss the consequence of imposing the condition of validity of the fully 
flavored regime, Eq. (36). The result is shown in the left panel of Fig. 9, to be compared 
with Fig. 6 where the same parameters were varied but all points were kept. One notices 
clearly that many points at large mi disappear and an upper bound on mi, given by 
mi < 0.15 eV, appears again and is comparable to the upper bound holding in the 
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mJeV] 



Figure 9: Same as Fig. 6, imposing the condition of validity of the fully flavored regime 
(left panel) and the strong washout condition (right panel). 



unfiavored regime, rrij < 0.12 eV. In the case of inverted hierarchy similar results apply. 
However, it should be clarified that mi < 0.15 eV is not an upper bound on the neutrino 
masses but a limit of validity of the fully flavored regime. One should therefore solve more 
general kinetic density matrix equations in order to describe the intermediate regime and 
to see whether there is or not an upper bound on neutrino masses. In the right panel 
of the same figure, we also imposed the strong washout regime condition and one can 
see that the allowed region gets further reduced and the upper bound on mi even more 
stringent, mi < 0.10 eV. 

So far all results have been obtained making use of the approximation C = I. In 
the upper panels of Fig. 10, left for normal hierarchy and right for inverted hierarchy, we 
show the results without taking into account the Higgs asymmetry solving the Eq.'s (24) 
neglecting C^ in C that reduces to C' given by the expression (25). One can see that the 
upper bound on mi gets relaxed when C^ is neglected. In the bottom panels we used the 
Eq. (27) for C and one can see that the results agree very well with those obtained when 
we used the approximation C = I. Therefore, we will continue to use this approximation 
in the following analysis. 

A particularly interesting case to be discussed is when all the high-energy phases are 
switched off, i.e. imposing the Q matrix to be real. In this case the PMNS phases act as 
the only source of CP violation responsible for the explanation of the observed matter- 
antimatter asymmetry [10]. It is therefore interesting to understand whether there is an 
allowed region at all [6]. The results are shown in Fig. 11 for M3 = M2 = 3 Mi. We show 
both the results obtained when both Majorana and the Dirac phase are switched on (red 
points) and those obtained when only the Dirac phase is non- vanishing (green points) and 
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Figure 10: Check of the approximation C = I. In the upper panels we neglected the 
Higgs asymmetry solving the Eq.'s (24) with C^ = 0. In the bottom panels we solved 
the Eq.'s (24) with C given by the Eq. (27). 
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for sin ^13 = 0.2, the 3 a experimental upper bound. This even more special case is the 
Dirac phase leptogenesis scenario [19]. The Dirac phase being the only phase that we can 
realistically measure in the near future in neutrino oscillation experiments, it is therefore 
particularly relevant to understand whether this testable source of CP violation can be 
also responsible for the observed matter-antimatter asymmetry. 

The top-left panel is for normal hierarchy while the top-right is for inverted hierarchy. 
One can notice that in this last case the lower bound on Mi tends to infinite for vanishing 
TTii, meaning that the final asymmetry tends to vanish ^. In both cases there is an upper 
bound on mi, as first discussed in [17]. In the bottom-left panel we impose the condition 
of validity of the fully flavored regime given by the Eq. (36) and one can see that the 
allowed region gets reduced in a relevant way. Imposing even the strong washout regime 
condition (i.e. independence of the initial conditions) one can see that the allowed regions 
gets further reduced and it practically disappears for the case of Dirac phase leptogenesis, 
confirming, more generally, what was found in [19]. 

5.3 Effect of large \uij\ values 

So far, we imposed the restriction \ujij\ < 1. In this subsection we relax this restriction 
studying the effect of allowing \ujij\ as large as 10 on the bounds. The main effect is 
that the extra-term Aeia in the flavored CP asymmetries (cf. Eq. (52)) can become 
dominant. The upper bound Eq. (55) does not apply to this term. Since this term is 
suppressed like M1/M2 the dominance is possible for a mild hierarchical spectrum, such 
that M2/M1 < (9(100). In any case notice that its dominance applies under conditions 
that are much less restrictive than those found for the dominance of A^i, that is suppressed 
hke (Mi/Ms)^ and that cancels when M2 = M3. 

In Fig. 12 we show the allowed region in the plane {Ki,Mi) for mi = and for 
M3 = M2 so that A^i does not give any contribution. The results have to be compared 
with those in Fig. 4. The maximum effect is obtained again, like for A^i, when 00^2 is large. 
In the top panels only |a;32| is allowed to be as large as 10, while |cc;2i|, |a;3i| < 1. In the 
top-left panel and in the top-center panel, for normal and inverted hierarchy respectively 
and for M2 = 3 Mi. One can see how the lower bound on Mi gets relaxed of one order of 
magnitude so that it can be as low as 10^ GeV. Analogous relaxation applies to the lower 
bound on the reheat temperature Tj-ch- Notice that in the limit of no washout, for Ki -^ 0, 
the usual unfiavored lower bound is recovered. In the top-right panel M2 = 30 Mi, such 
that the extra-term Aeia is suppressed enough not to be able to relax the lowest bound 



^Similar results have been recently obtained in [47]. 
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Figure 11: Leptogenesis only with Majorana phases (red points) and only with Dirac 
phase (green points). Top- left and top right panel: allowed region in the plane (mi, Mi) 
for normal hierarchy and inverted hierarchy respectively. Bottom panels: the same as 
top-left but imposing the condition Eq. (36) (left) and the strong washout condition as 
well (right panel). 
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Figure 12: Lower bound on Mi versus Ki for rrii = 0. The results are compared when 
flavor effects are neglected (green points) and when they are taken into account (red 
points) while all have been obtained for a thermal initial A'^i-abundance (A?^]^^ = !)• Top 
panels: |ci;2i|, |ct;3i| < 1, |u;32| < 10; Top left and center panels: normal and inverted 
hierarchy, respectively, and M2 = M3 = 3 Mi; Top right panel: normal hierarchy and 
M2 = M3 = 30Mi. Bottom panels: {uj^il, \uj32\ < 1, |i^2i| ^ 10; Bottom left and center 
panels: normal and inverted hierarchy, respectively, and M2 = M3 = 3 Mi; Bottom right 
panel: inverted hierarchy and M2 = M3 = 30 Mi. 
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on Ml. In the bottom panels we allowed |c<j2i| to be as large as 10 while \uj3i\, \uj32\ ^ 1- 
One can see that the effect is much more reduced since now A^iq, increases but at the 
same time Ki and the washout increase as well. 

6 On the validity of the A^i-dominated scenario 

In general, the final asymmetry receives a contribution from the decays of all 3 RH 
neutrinos: N^b-l ~ Xli ^b-l(^«)- -'■^ ^^^^ Section we want to find a condition that 
guarantees the validity of the bounds within the iVi-dominated scenario, where N^b-l — 

N^B^LiNl). 

If M3 > 10^^ GeV, then the asymmetry produced from the heaviest RH neutrino 
decays occurs always in the unfiavored regime. If moreover we impose M3 ^ M2, then it 
can be always safely neglected since the total CP asymmetry e^ oc {Mi^2/M^)'^ is strongly 
suppressed and at the same time a strong washout from A^i and N2 is unavoidable [5] . 

On the other hand, 62 is not necessarily suppressed. A condition that guarantees 
the possibility to neglect N^^_j^{N2) in the determination of the bounds is equivalent to 
impose 

r/B(iV2)^a.p,^%^«^r^ (58) 

From the Eq. (8) and using the orthogonal parametrization, the CP asymmetry 62 can be 
recast as [5] 



62 = e{M2 
where we defined 



4Mf /M2 
3 M^ \Mi 



(59) 



/923(mi,fi) + - ^ In ( ^ - 1 ) P2i{m,,n) 

A,K,(])^ME^!!^^l!Mm1!. (60) 

rrii matm 

There are two special cases for which the A'^i-dominated scenario is certainly valid. The 

first is to have Q = R13 (i.e. a;2i = 1^32 = 0), since in this case it is simple to see that 

82 = 0. The second case is the limit M3 ^ 10^^ GeV. In this case Q is given by the 

Eq. (33) and it is easy to see that £2 = 0. This result can be understood considering that 

in this limit the heaviest RH neutrino decouples and so necessarily the interference term 

(m|)mz))23 -^ 0. 

The opposite case is for ^2 = i?23, since now one has £1 = while 

82 < e{M2) ""''""' /(m2, m2) , (61) 
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the same maximum value holding for ei (cf. Eq. (49)) but where (mi, Mi, mi) are replaced 
by {1712, M2,rh2). Notice moreover that for Q = i?23 one has mi ^ m^, so that the 
asymmetry produced from A'"2-decays is certainly not washed out by A''i-inverse decays. 
In this situation a A''2-dominated scenario is realized, with N^^-l — ^B-ii-^^) ^^^ with 
no lower bound on Mi [5]. 

Between these two well-defined special cases, one has to take into account both a 
contribution iV^_^(A^i) from A^i-decays and a contribution A''^_^(A^2) from A'^2-decays. 
For example, choosing Q = R12 one has 

4 Ml 

Assuming that M2 ^ 5 x 10""^^ GeV, such that the asymmetry from A^2-decays is pro- 
duced in the unflavored regime and that M2 ^ Mi, one can see that even neglecting 
the washout from A''i-inverse processes the contribution N^q_j^{N2) from A'"2-decays can be 
safely neglected in the determination of the bounds. 

This example shows that if the first term in the Eq. (59) oc lra[{m\j mi)) 23]"^ vanishes, 
then riB{N2) can be neglected in the determination of the bounds in the hierarchical 
limit where M2 ^ Mi. Therefore, a condition Im[(m})m£))23]^ = certainly guarantees 
the validity of the A'^i-dominated scenario but is quite a restrictive one. However, this 
condition enlightens that an asymmetry generated from A^2-decays requires an interference 
of the heaviest RH neutrino A^3 with A^2- Indeed, as we said, this condition is certainly 
verified in the limit M3 ^ 10^^ GeV, when the heaviest RH neutrino decouples and 
(m^mz))23 = 0. 

Now we want to see whether, starting from Vt = R13, one can turn on a rotation 
-R12 (^21 7^ 0) still having a negligible rjB{N2). Since both for Vt = R13 and Q = R12 
one has {m.\jm.£))23 = 0, one could naively think that this is still true for n = R12R13 
and therefore that a;32 = is a sufficient condition for the validity of the A'^i-dominated 
scenario. However it is easy to check it is not true that {m}j^mbj:,)2'i = and therefore 
it is not guaranteed that the asymmetry from A^2-decays can be neglected. This has to 
be done by inspection. Since we are assuming the hierarchical limit, M2 > 3 Mi, the 
calculation of N^^-ii^'i) factorizes in two terms 

Vb{N2) = r]B{N2)\T^T^^K,) X MT ~ Ml) . (63) 

The first term is the asymmetry produced at Tg ~ M2/ zb{K2) from A^2-decays, while the 
second term is the washout from A'^i-inverse processes. 

A precise calculation of Wi{T ~ Mi) has to take into account two types of flavor 
effects. If Ml < 10^ GeV, then the asymmetry produced from A^2-decays is fully projected 
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on the three-flavor basis at the time when the asymmetry is washed out by A'^i-inverse 
processes [49]. Assuming moreover that M2 ^ 5 x lO^"*^ GeV, so that the asymmetry from 
N2 decays is produced in the unflavored regime, the contribution to the final asymmetry 
from iV2-decays can be calculated as 

NUiN^) = e, k{K,) J2 PL e-t <^^ . (64) 

a 

On the other hand, if M2 ^ 5 x 10^^ GeV, then the asymmetry is produced in the fully 
flavored regime and 

Ni_^{N,) = 5^e2««:(i^2«)e-^<^^ . (65) 

a 

Notice, however, that since we are interested in finding the condition for the A'^i-dominated 
scenario to hold, then Mi > 10^ GeV and at the time of the washout from A'^i-inverse 
processes the asymmetry is projected on a two-fiavor basis: the r fiavor and an orthogonal 
combination of electron and muon fiavors. In this situation A^i-inverse processes have the 
effect to further project the asymmetry stored in the e+/i flavor on a two-flavor basis where 
one direction is determined by /i and the other is the orthogonal component [12, 48]. The 
washout from A^i-inverse decays does not touch this orthogonal component and therefore, 
accounting for this effect, the value of the flnal asymmetry has to lie somewhere between 
the value of the asymmetry produced at T ~ Tb{K2) and the value of the asymmetry 
calculated neglecting the effect of projection along £1, 

NB-dN2%^^^^^,) X wiiT ~ MO < N'B_dN,) < Ni_,{N,)\^^^^^^^^ . (66) 

As already mentioned, two special cases where the validity of the A'^i-dominated scenario 
is guaranteed are M3 ^ M2 > 10^'^ GeV and Q = Ri^. In Fig. 13 we show the results 
obtained in a more general case. We allow Q = i?i2(^2i) -^13(^31) where 0^21 7^ but 
still 00^2 = 0. In the left panels we show the results in the plane {M2,i]b{N2)). The 
flnal asymmetry is calculated both neglecting the washout (red points) and with the 
washout calculated in the two-flavor regime (green points) but neglecting the effect of 
projection of the asymmetry due to A'^i-inverse processes: as we said an account of this 
effect should give a result that has to be somehow in between. One can see in the top-left 
panel that for M2 < 10^^^^^ GeV the asymmetry produced by A^2-decays falls below the 
observed one indicated by the horizontal solid line. Here we are imposing a;2i < 1 and we 
are assuming the reheat temperature to be higher than ~ Tb(A^2) ~ ^2/ zb{K2). This 
result can be actually also translated into a condition on the reheat temperature given by 
^rch ^ Tb{N2) ~ M2/ zb{K2)i while M2 this time is free. The result is shown in the right 
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MjGeV] 





Figure 13: final asymmetry from A''2-decays versus M2 (left) and versus the reheat tem- 
perature Treh (right) for |u;32| = and tusi < 1. In the top panels UJ21 < 1, while in the 
bottom panels UJ21 < 0.1. The final asymmetry is obtained neglecting the washout term 
Wi{T ~ Ml) in the Eq. (63) (red points), calculating it taking into account flavor effects 
(cf. Eq. (64)) (green points) and neglecting them (blue points). 
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panels of the same figure and one can see that the lower bound on the reheat temperature 
is approximately five time more relaxed compared to the lower bound on M2. 

In the bottom panels we imposed \uj2i\ ^ 0.1. One can see that even for such small 
angles there is still a marginal allowed region. Therefore, only for very small \uj2i\ values 
one recovers the special case Q = R13. 

This result shows that while for Q = R13 or for Q = R12 the asymmetry from N2- 
decays is always suppressed and the A'^i-dominated scenario holds, as soon as one allows 
either u^i 7^ or CJ21 7^ or both, the asymmetry produced from A^2-decays can explain 
the observed asymmetry for acceptable values of M2. Even though we did not perform 
a systematic calculation of the final asymmetry from A^2-decays in the whole parameter 
space, in the light of this result it emerges that, if 10^^ GeV < M2 < 10^^ GeV, then 
the iV2-dominated scenario seems to be the most natural choice and the validity of the 
iVi-dominated scenario relies on unnaturally small complex angles in order to suppress 
Im[(m}^ 771^:1)23]^. We can therefore confirm even in a stronger way what we said in the 
introduction: it is misleading to talk of a lower bound on Mi in leptogenesis, while it 
is more correct to talk of a lower bound on T^eh holding in the hierarchical limit for 
M2>3Mi. 

We also notice that this result is not relying crucially on an exact calculation of the 
washout from A'^i-inverse processes but more on the calculation of the asymmetry produced 
at T ~ M2. In the same Fig. 13 we show (red points) the asymmetry ^b(A^2)It~t (k )' 
without the washout, while the green points take into account the washout term Wi{T ~ 
Ml) calculated with the Eq. (64). We also show (blue points) the results that would have 
been obtained calculating the washout neglecting flavor effects. 

One can see that the main result relies on the observation that ^b(A^2)Ii-^t (k-) ^^ 
large already when small \uj2i\ are turned on. A proper calculation of the washout is an 
important step but somehow a secondary one. Certainly a proper account of flavor effects 
greatly enhances the asymmetry from A^2-decays but primarily it is important that there 
is an asymmetry to be washed out. Notice that an account of the effect envisaged in [48] 
would produce results somewhere in between the green and the red points. It is certainly 
important for a precise evaluation but it does not seem to play a crucial role. 

Our results show that there is a continuous increase of the asymmetry going from 
Q = i?i3, where it vanishes, toward Q = R23. What is important is that as soon as one 
switches on small \uj2i\ or {uj^il the contribution from A^2-decays is sufficient to explain the 
observed asymmetry. Going toward the case Q = R23 (u;2i = ^31 = 0), the allowed region 
for the iV2-dominated scenario increases and the lower bound on M2 relaxes down to a 
few's xlO^° GeV [5]. We do not see any discontinuity or qualitatively different regime as 
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envisaged in [22] where the authors distinguish a decoupled regime from a strong washout 
regime. A proper evaluation of the washout term Wi{T ~ Mi), taking into account 
all kinds of flavor effects [49, 48] is in any case an important ingredient for a correct 
determination of the border between the domain of validity of the A^i-dominated scenario 
and that one of the A^2-dominated scenario. 

7 Beyond the hierarchical hmit 

Finally, in this section we discuss how the bounds change going beyond the hierarchical 
limit, when 62 = (M2 — Mi)/Mi < 2. We first neglect flavor effects induced by charged 
lepton interactions. Notice that there is a second type of flavor effects in the heavy 
neutrino sector itself [12, 24, 48] due to the fact that for example the lepton quantum 
state 1/2) produced by a RH neutrino N2 does not coincide with |/i) produced by a RH 
neutrino A^^i and in particular I/2) does not in general inverse decays with a Higgs to 
produce a A^i with the same rate as |/i). We will neglect this effect that would imply that 
washout terms do not simply add up in the Eq. (16). 

Under these assumptions, there are three different effects that change the bounds 
compared to the hierarchical limit [43]. First of all, in general, one cannot neglect the 
contribution from the heavier RH neutrinos. Therefore, one cannot just relax the assump- 
tion of hierarchical spectrum without also considering the heavier RH neutrino decays. 
Second, now the addition of washout in the Eq. (16) cannot be neglected. Third, the 
total CP asymmetries get typically enhanced compared to their value in the hierarchi- 
cal limit. The first and third effect tend to increase the final asymmetry relaxing the 
bounds but the second effect tends to reduce the final asymmetry, making the bounds 
more stringent. The first two effects saturate for 62 ^ 0.01, the so called degenerate limit, 
and therefore when 62 decreases below 0.01 only the third effect is left and it changes the 
bounds in quite a simple way [43], except for some effects studied in [50] in the extreme 
case of resonant leptogenesis. However, notice that in the case of resonant leptogenesis the 
bounds simply disappear. Therefore, here we focus especially on the transition between 
the hierarchical limit and the degenerate limit, for 0.01 ^ ^2 ^ 2. The final asymmetry 
is given by N^b~l ~ Xlj ^« '^i) where ei is given by the Eq.'s (37)-(40) while, from the Eq. 
(8), the total CP asymmetries £2 and £3 are given by 

^ 3 ^p Im [{h^h)%] i{x,/x2) ^ 3 ^ Im [{h^ h)%] i{x,/x^) 

^2 16^ Z. (;,t/,)22 ^/^ ^"^ '' 167r^ {WhU V^' 

(67) 

35 



The unflavored efficiency factors are given by the Eq. (35) with P^^ = 1 



where we defined 6^ = (M3 — Mi) /Mi. Let us assume that the A^j-abundances track 
closely the equilibrium value, so that dN^Jdz ~ dN'^Jdz. This approximation works 
well for fTj > 1, as we will assume in the following. An approximated expression for kI is 
then obtained 



kI{Kj,62,63] 



dz' -^ expl- / dz" [A Wiz") + J2 ^n^")] \ ■ (69) 



In [5, 43] it was shown that conservatively for 

^i ^ ^HL = —j-^^ — ^ - 1 , (70) 

Zb{j<i) -^ 

one can assume the hierarchical limit for the calculation of nj. In this case the washout 
from the RH neutrinos / 7^ i lighter than Ni is factorized and (cf. Eq. (44)) 

^l ^ n{K,) e- /o°° '^^' S' ^"(^') . (71) 

On the other hand the washout from the iVj-inverse processes on the asymmetry produced 
by iV^-decays is negligible and 

^l ^ _ r dz' ^^ e- f^ '^" [^^(^'O+E, w}-i,.";m , (72) 

where with A*";^ we indicated one particular A";. Typically one has (5hl — 2 and this is the 
reason why we always used M2 > 3 Mi as a condition for the hierarchical limit to hold. 
On the contrary, if (Mj — Mi)/Mi < 0.01, then one recovers the degenerate limit where 
kI ^ f^\{Ki + Ki,KjjLii). There are three different cases: a partial degenerate limit with 
62 ^ 0.01, a partial degenerate limit with (M3 — M2)/M2 < 0.01 and a full degenerate 
limit where ^3 < 0.01. In this last case one has simply 

4 = 4 = «3 - '^(^1 + K2 + Ks) . (73) 

We will now focus on two particular choices for M3 but without imposing any restriction 
on 62- The first case we consider is M3 ^ M2 ~ Mi. This was also studied in [43] 
but for two particular choices of the orthogonal matrix, while here we do not make any 
assumption on Q. The contribution from the heaviest RH neutrino is negligible both 
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because /t| <^ k{, k\ and because £3 <^ ^i, £2- The washout of A^3-inverse processes on the 
asymmetry produced from the two hghtest RH neutrino decays is neghgible as welL 
Two convenient fits for k\ and k\ can be used [43], 

kI\Ki, K2, 62) = ^^4 r- (74) 



and 



M 



I-S2 ^ 



(K.K,M = .}!!-:''. ..._.,, X e--'MA) . (-5) 



For M3 ^ M2 the Eq. (67) for 62 can be further speciahzed into 

m2 matm "^2 ^atm V 2^2 / 

If ^2 ^ 1 and if one maximizes over the Q parameters, the first term can be neglected 
and since .^(l/x2) — —^(2^2), the expression simphfies into 



maxf7[£2] ^ -£{Mi) ^(^s) /92("^i) , P2{mi) = max^ 



Im [J2h rrih ^1,2^ 



hl\ 



l2 



^2 matm 



(77) 



While El gets suppressed when mi increases (cf. Eq. (49)), the function j32{'mi) in general 
does not and therefore, at large mi, it gives the dominant contribution to the maximum 
asymmetry determining the bounds in the plane (mi. Mi) that are shown in Fig. 14 for 
62 = 1,0.1,0.05,0.01 and where we imposed a strong washout condition. This automati- 
cally guarantees the validity of the fits (74) and (75) for the efficiency factors. The results 
in the unfiavored case correspond to the green points. One can notice a result already 
found in [43] that we confirm here in a more general way. For ^2 ~ 0.1 the washout addi- 
tion makes the lower bound on Mi even slightly more stringent while only for 62 ^ 0.01 
the lower bound gets clearly more relaxed. On the other hand, concerning the upper 
bound on mi, the fact that 62 increases with mi implies that the upper bound on mi gets 
relaxed already at 62 — 0.1. 

Let us now turn to study the second case when M3 = M2. This time the heaviest RH 
neutrinos contribute both to the asymmetry production and to the washout. Having the 
expressions for the case M3 ^ M2, it is easy to derive the efficiency factors in the partial 
degenerate hmit (M3 — M2)/M2 < 0.01. Indeed if i = 2, 3, one has simply 

4^4 ^ -y d^ ^ X (78) 

X exp |- r dz" [A W{z") + W^^{z"- Ki) + WJ^'iz"- K2 + ^^3)] 

~ kI\Ki,K2 + K^,62). 
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Figure 14: Allowed region in the (mi, Mi) plane for the case M3 ^ M2 in the unfiavored 
case (green points) and in the fully flavored regime imposing the condition Eq. (36) (red 
points) for ^2 = 1 (top-left), 0.1 (top-right), 0.05 (bottom-left) and 0.01 (bottom-right). 
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The calculation of the total CP asymmetries £2 and £3 is slightly more involved. If we 
assume that M2 = M3 exactly, then the interference term between N2 and N^ vanish, 
though notice that the expression Eq. (8) for 82 and £3 diverge for M2 = M3 because they 
hold only for mass differences less than the decay widths. Therefore, one has 



ei=2,3-^{x2)€{M,) 



E. ml lm[f]*,2 ^2 J ^^y^^^^^ Im[fiM Qn ^l, ^t 



T^i TTlatm 



h<l 



mii Tnatm 



(79) 



In passing let us notice that, without the interference term between N2 and A^3, the 
asymmetries vanish for Mi -^ 0. The final asymmetry can then be written as 



N- 



B-L 



e(x2) 5(Mi) {nf{K,, K2; 62) (3{m,, Q) + 4\K,, K2 + K^, 62 



iO) 



X 



2 n 



j:H<M^f2^li] , En^M^fs^li 



+2 ^ mhmi 

h<l 



m2 rriatm m^ m^tm 

Im[Qhi ^11 ^Is ^h] ffh2-m^ 



^atm 



1712 1713 



(81) 



Notice that if m2 = "^3, then £2 + ^3 = (^2/'^i)^i oc {m^ — rrii) and therefore the 
contribution from the two heavier RH neutrinos is suppressed as well. This makes the 
bounds more stringent compared to the case M3 ^ M2, especially for 62 ~ 0.1. In Fig. 15 
we show the bounds in the case M2 = M3 again for 62 = 1,0.1,0.05,0.01. One can see 
that now for ^2 ~ 0.1 not only the lower bound on Mi does not get relaxed but also the 
upper bound on mi (green points). Again, for 62 ^ 0.01, both the Mi lower bound and 
the mi upper bound get clearly relaxed. 

Both in the case M3 ^ M2 and in the case M3 = M2, we repeated the calculations 
also in the fully fiavored regime. Now the final asymmetry has to be calculated using 



N 



B-L 



E 






Y,N'f,e,Kl + lY.AP,^[Kl 



^iidl ■> 



12) 



where we generalized the approximated form we already used in the hierarchical limit (cf. 
Eq. (50)). 

In the case M3 ^ M2, the contribution from A''3-decays is this time less trivially 
negligible because in principle the washout from the the two lighter RH neutrinos could 
be weaker in one fiavor and because moreover the fiavored CP asymmetries £30,, contrarily 
to the total £3, are not necessarily suppressed in the hierarchical limit [6]. However, it 
turns out that the contribution from A'3-decays is too small to explain the observed 
asymmetry and it can be therefore neglected in the determination of the bounds. 
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Figure 15: Allowed region in the {mi, Mi) plane for M2 = M3 in the unflavored case (green 
points) and in the fully flavored regime imposing the condition Eq. (36) (red points) for 
62 = 1 (top-left), 0.1 (top-right), 0.05 (bottom-left) and 0.01 (bottom- right). We are also 
imposing \u!ij\ < 1 and K^i > 1. 
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Assuming Kai ^ 1, expressions for the flavored efficiency factors can be obtained 
from the unfiavored ones just replacing Ki — > Ki^ = P^^Ki. In this way we can use the 
approximation (69) and the fits Eq. (74) and (75) with Ki replaced by Kia. 

In Fig. 14 one can see how flavor effects affect the bounds (red points). For ^2 = 1 
there is no relaxation of the lower bound on Mi, since the hierarchical limit still holds and 
because at the onset of the strong washout there is no relaxation due to flavor effects for 
\ujij\ < 1 [6], as we are imposing. For ^2 ^ 1 and mi < 0.01 eV one can see that there is a 
factor two relaxation. The reason is that the asymmetry is not maximized in a one flavor 
dominance case but in a democratic case, that means for values of the parameters where 
the only change to the final asymmetry from flavor effects is described by the enhancement 
of a factor N"^ ~ 2, while the additional terms oc APjq, (kjo, — nip) vanish. On the other 
hand, for mi > 0.01 eV, one flavor dominance makes possible a large enhancement of the 
asymmetry and this is why the upper bound on mi is much more relaxed. Similar results 
hold in the case M2 = M3, as one can see from Fig. 15. 

8 Conclusions 

The simple vanilla leptogenesis scenario grasps important features of leptogenesis bounds 
but misses many important effects. Assuming the A^i-dominated scenario, the lower bound 
on Ml seems to be a solid feature and we have seen that it resists even for heavy neutrino 
mass degeneracies as small as 82 ~ 0.01. However, our analysis revealed that flavor 
effects introduce new CP violating terms that relax the bound of one order of magnitude 
for acceptable choices of the parameters and still within the hierarchical limit. Flavor 
effects modify the upper bound on mi as well but, as we stressed, an ultimate answer 
requires solutions of more general kinetic equations that should be able to describe the 
intermediate regime where the coherence of the final quantum lepton state is lost but a full 
decoherence is still not achieved. We have seen that an account of the Higgs asymmetry 
supports such a prudent conclusion. 

Still within the hierarchical limit but accounting for the asymmetry produced from 
the A^2-decays, the lower bound on Mi disappears and is replaced by a lower bound on M2 
that still implies a lower bound on the reheat temperature. We showed some more general 
conditions for the validity of the A^i-dominated scenario and of the lower bound on Mi. 
For example the iVi-dominated scenario certainly applies in the popular two effective RH 
neutrino limit, for M3 ^ 10^^ GeV, that implies that the heaviest RH neutrino decouples. 
However, apart from this case, our analysis indicates that a A^2-dominated scenario is a 
more natural option and that neglecting the asymmetry from A^2-decays can be a wrong 
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assumption. 

In conclusion, our analysis answered many different questions about the bounds on 
neutrino masses that are obtained within the leptogenesis scenario based on the simplest 
version of the see-saw mechanism. Despite many proposed extensions, this still represent 
the most attractive possibility since it realizes a successful link between the neutrino 
masses and the observed asymmetry where the measured values exhibit an interesting 
conspiracy. The discovery of CP violating effects in neutrino oscillations or in lepton 
decays, the determination of the absolute neutrino mass scale and of the neutrino mass 
spectrum ordering, normal or inverted, will likely provide further interesting tests during 
next years, making current experimental 'coincidences' even stronger or forcing departures 
from the minimal picture. 
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